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Abstract 

We discuss the use of the variational principle within quaternionic quantum mech- 
anics. This is non-trivial because of the non commutative nature of quaternions. 
We derive the Dirac Lagrangian density corresponding to the two-component Dirac 
equation. This Lagrangian is complex projected as anticipated in previous articles 
and this feature is necessary even for a classical real Lagrangian. 
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1 The variational principle 



The standard derivation of the Euler-Lagrange field equations implicitly assumes that 
the variation in the action / of a given field Sfa may be commuted to, say, the extreme 
right. Thus after a functional integration by parts, and neglecting the surface terms 
one obtains the field equations, 



where C is the Lagrangian density and x represents space-time. 

Somewhat surprisingly this form of the field equations survives in field theory, 
that is after second quantization. The reason that this is not a priori obvious is 
because when the fields become operators and non-commutating, the translation of 
5 fa is not necessarily without consequences. Indeed, we should be obliged to define 
with attention the significance of the functional derivatives in eq. (|l|). This fact does 
not in practice produce difficulties as a consequence of: 

1. For bilinear terms (e.g. mass or kinetic energy) the variations naturally lie on 
the right (for 5 (pi) or the left (for Scf)f); 

2. In any case some authors implicitly assume that, when 5 fa is bosonic it commu- 
tes with everything; 

3. When fa is a fermion field 5<pi is assumed at least to commute with any bosonic 
fields and with fa itself. This last fact (by no means obvious) leads to a natural 
extraction of 5 fa to the right or the left. 

More specifically, within normal ordered products there always exists a natural dire- 
ction of extract given the fact that all creation operators ( bosonic or fermionic) 
commute amongst themselves. 

It is however an assumption that the variation in fa, i.e. Sfa has the same char- 
acteristic commutation relations of fa. These implicit assumptions are brought to 
the fore when one considers quaternionic fields both in classical and quantum field 
theory. For it is unthinkable that a quaternionic 5<f)i created within a Lagrangian 
density can, in general, be commuted without consequence to the left or the right. 
The example of the free Dirac Lagrangian treated in this paper will demonstrate some 
of the difficulties. 

Before entering explicitly into the world of quaternionic fields we wish to discuss 
the limitations, if any, applied to 8 fa in standard (complex) quantum mechanics. Our 
objective, which will have an application in this paper, is to demonstrate that the 
properties of 8 fa may be substantially divers from those of fa. For example, we shall 
observe first that even if the classical Lagrangian density is necessarily real, 8<pi or 
more correctly 8C may be formally complex (or even quaternionic) since at the end 



dC{x) 
dfa(x) 



0, 



dCjx) 
dfj,fa(x) 



(1) 
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5C = to obtain the Euler-Lagrange equations. It is true that in general we may limit 
our attention to a subset of 5(f>i (S£) without losing the Euler-Lagrange equations, 
so that our observation seems accademic, but it is not without consequences with 
quaternions as we shall see. 

Consider one of the simplest of all particle Lagrangian densities, that for two 
classical scalar fields (0j-real i = 1, 2) without interactions: 

c = - — # + ~w> 2 - — 2 2 

= d^d^cj) - m 2 (2) 



where 



and 



+ k 



l(p2 



The well known corresponding Euler-Lagrange equations are: 

(d^ + m 2 )0 i = O , (3) 
(* = 1, 2) , 

or, equivalently, 

(d^ + m 2 )(j) = . (4) 

Now to obtain "directly" the last equation one performs the very particular variation 
of 0(0+) 



0+ -> (p + + 5(j) + (5) 

i.e. in order to obtain the corresponding Euler-Lagrangian equation one treats (j) and 
(f) + as independent fields. In second quantization these fields indeed contain inde- 
pendent creation and annihilation operators corresponding to positive and negatively 
charged particles^. To satisfy eq. (|5|) we must necessarily have, 

Sfa + z50 2 = (6) 

and this means, that the variations in the originally real (pi fields are complex (if 5(fii 
is real then 5(p2 is pure imaginary etc.). 



1 Actually the "electric" charge e of the fields depends upon whether the global gauge group is to 
be made local or not. This fact is a choice not determined by the free Lagrangian. 
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Of course we could obtain the equivalent result from variations separately of 0i 
and 02, e.g. varying 0i (with 2 constant) 

0^0 + 501 
0+ - 0+ + 50! (7) 

yielding, after a double integration by parts 

(5 M ^ + m 2 )(0 + 0+) = (8) 

and the corresponding result for 2 . Thus while not obligatory we can in this latter 
approach consider only real <50j. However, there is a subtle difference in the two 
approaches, which readily passes unobserved. In the case when only 0i or 02 is 
varied the 5<pi appears both to the left and to the right. Only if 50, commutes with 
4>i are the two methods equivalent. 

The point of these observations is that it is possible to impose diverse conditions 
on the Lagrangian density, on the fields, and on their variations. This is ignored in 
standard quantum mechanics but is important for what follows. 



2 The two component Dirac Equation 

The adoption of quaternions as the base number system in quantum mechanics allows 
one to define a quaternionic Dirac equation |]J in which the wave function is charac- 
terised by having only two components (consequently the Schrodinger-Pauli equation 
applies to a one-component wave function). This follows from the fact that the Dirac 
algebra upon the reals (but not upon the complex) has a two dimensional irreducible 
representation with quaternions. The standard 4x4 complex gamma matrices are in 
fact reducible. This structure is consistent if the momentum operator p^ is defined 
by 

= | % (9) 
where a barred operator A \ b acts upon a general wave function ip by 

(A | b)ifj = Ai\)b . (10) 

In general A will be a matrix and b a complex C(l, i) number, where i is one of the 
imaginary (i, j, k) units of a quaternion: 



with 



q = a + aii + a 2 j + a 3 k , 

i 2 =f = e = -l , 
ij = —ji = k ( cyclic ) , 



(11) 
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Oq, 1, 2, 3 £ 1Z . 

The quaternionic conjugation is denned by 

q + = cio — a\i — a 2 j — a 3 k . (12) 
The Dirac equation then reads in covariant form 

j^tpi -rmf; = (13) 
with as a possible choice of the 7 M 

70= (J -l) ' 7l = *(? o) ' 

■>>-> ■(:;).*- (i4) 

This equation has four plane wave (oc e _4p:r ) solutions which correspond to those 
of the standard Dirac equation if one adopts the complex scalar product (complex 
geometry ||) 

1 % 

<i/)\(l)> c ^-<ip\(j)>--<ip\(j)>i (15) 

where < ip \ 4> > is the quaternionic scalar product, 

<if;\<P>= J ifj + (f) dr (16) 

so that the four plane wave solutions are (complex) orthogonal. Actually the need 
of this scalar product is anticipated by our choice of which is not hermitian on a 
quaternionic geometry. The complex scalar product was first introduced in the defini- 
tion of tensor products in quaternionic quantum mechanics 0, and appears essential 
for the "translation" of complex quantum mechanics to a quaternionic version [Q. 

Our main objective in this work is to derive the Dirac Lagrangian (density) which 
yields eq. flT3]), and this will be done in the next section. We conclude this section 
with some notes and comments upon this use of quaternions in quantum mechanics. 

1. The objective of reproducing almost all standard results in quantum mechanics 
is achieved; 

2. There are differences within the bosonic sector. For example there is a doubling 
of solutions in the Klein-Gordon equation 0; 

3. There exist however projected equations with only the standard number of 
solutions [Bl |7f|; 
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4. A particular example of point 3 is the non equivalence (with quaternions) of the 
Maxwell equation with the Dumn-Kemmer-Petiau equation. The latter has no 
doubling of solutions, and may be identified with a projected Maxwell equation; 

5. It is possible to invent new quaternionic equations equivalent to pairs of complex 
equations in the same way as the Schrodinger equation can be rewritten as a 
pair of real equations; 

6. There exists a not yet completely investigated generalization of group theory or 
more precisely of representation theory. 

Finally, we are always intrigued by the fact that had Schrodinger considered 
quaternionic solutions to his equation, he would have found two (within a complex 
geometry) and probably discovered the existence of spin. 



3 The Dirac Lagrangian 

Let us consider, as a first hypothesis, the traditional form for the Dirac Lagrangian 
density: 

C = l -[4>Y{d^) - M-fif*] - (17) 

as given for example by Itzykson and Zuber ||. The position of the imaginary unit is 
purely conventional in ( |TTD but with a quaternionic number field we must recognize 
that the <9 M operator is more precisely part of the first quantized momentum operator 
<9 M | % and that hence only the hermitian conjugate part of the kinetic term (second 
half of the bracket expression in eq. (|T7|)) appears with the imaginary unit in the 
"correct" position. Thus the correct form of the kinetic term reads: 

C K = ]pYd^i-i{d^) 1 ^} . (18) 

We observe that this modification of eq. ( |TTD is also justified by the simple requirement 
that C be hermitian. We could of course multiply eq. (|18D by the "hermitian" 
which inverts the order of the imaginary unit, but by integrating by parts this may 
be reformulated as in eq. flI8|). Notice that we cannot change the sign of C without 
changing the sign of the Hamiltonian 7i. 

The requirement of hermiticity however says nothing about the Dirac mass term 
in eq. (17|). It is here that appeal to the variational principle must be made. A 



variation 5tf> in if) cannot in eq. ( |T8| ) be brought to the extreme right because of the 
imaginary unit in the first half of the expression. The only consistent procedure is to 
generalize the variational rule that says that if) and if) must be varied independently. 
We thus apply independent variations to if) (Sip) and if>i (5(if)i)). Similarly for Stf) and 
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5 (ill)). Now to obtain the desired Dirac equation for ip and its adjoint equation for ip 
we are obliged to modify the mass term into 

Tfl — — 

C m = - — [ipip - iipipi] . (19) 

The final result for C is 

I — — Tfl — — 

Cd = -{^d^i - i(d^)j^] --[W- iipipi] . (20) 

Considering this last equation we observe that it is nothing other than the complex 
projection of equation ([17]) 

Co = 1 ~^C . (21) 

Indeed, while C in eq. (|17D is quaternionic and with the modification of Ck in eq. ( |18D 
hermitian, the form given in eq. ( p0| ) is purely complex and hermitian. 

This observation is even more subtle with classical fields for now C defined by 

C = C K - m^ip (22) 

is both hermitian and real. Thus it may be objected that the complex projection 
in eq. fl2lD is superfluous. For C itself this is true but for quaternionic variations 
in the fields dip etc. a difference exists. The variation 5C of eq. (|2"2] ) is in general 
quaternionic while, because of eq. (|21|), the variation SCd is necessarily always com- 
plex. Furthermore eq. (|22| ) would not yield the correct field equation through the 
variational principle unless we limit dip to complex variations notwithstanding the 
quaternionic nature of the fields. We consider this latter option unjustified and thus 
select for the formal structure of £ that of eq. ( |20| ) . 

We note that as described in the first section, one should as a rule define the 
numerical basis of the fields and separately of the variations. For classical fields we 
have the case of real fields with complex variations (charged scalar fields) and in our 
studies the possibility of quaternionic fields but with only complex variations (the 
second option described above). 

Furthermore one should specify the numerical nature of both £ and 5C For 
example with our choice given by eq. (|2~0"[ ) we have: 

if) , tp , quaternionic fields ; 

Sip , 6(4>i) , Sip , 6(iip) , quaternionc variations ; 

Cd , hermitian, complex projected Lagrangian ; 

SCd , complex but generally non hermitian variation . 

4 The invariance groups of Cd 

Having obtained Cd in the previous section we may ask which if any group (global) 
leaves this Lagrangian invariant. Remembering that the single particle fields are 
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one-component quaternionic functions, the most natural transformation is, 



if) -> fipg 

$ - 9 + ^ + f + l° ■ (23) 

Now the quaternionic nature of the 7 M matrices limits /, and for arbitrary if) (since 
we assume no knowledge at this stage of the field equations^) leads to the conclusion 
that / must be real and commute with each 7^. On the other hand the complex 
projection of Cd permits g to be complex (it can then be "extracted" to the right 
and then commuted through Cd). Furthermore for invariance, 

99 + ff + = \9?f = 1 ■ (24) 

We find with our representation for 7^ that / is proportional to the unit matrix, and 
thus being real its magnitude may be completely absorbed within g, or equivalently 
/ may set by definition to 1. Whence, the only invariance group is defined by 

\9? = 1 , 

9 e C(l, i) . (25) 

For Lie groups this implies that g is a member of U(l, c) the complex unitary group. 

Remembering that the Glashow group || for the Weinberg-Salam theory [10| is 
SU(2, c) x U(l, c) we observe that this U(l, c) group may be identified with ours 
and that our field must necessarily be a singlet (scalar) under SU(2, c). Since if) 
represents a single fermion spin 5 field (the 7^ are 2x2 quaternionic matrices) it is 
not surprising, and indeed necessary that it represents a singlet under SU(2, c). 

The interesting feature is what happens if we select a reducible 4x4 complex set 
of 7 M matrices. Now the number of the fermion particles is two. Nevertheless, the 
existence of a complex SU(2, c) group acting on if) from the left seems excluded. On 
the other hand, as we have described in detail elsewhere |L1| the group SU(2, c) is 
isomorphic (at the generator algebra level) with U(l, q), the unitary quaternionic 
group with elements 

g ~ e ia+jb+kc 

with 

a, b, c G TZ , 

which is the simplest of all unitary quaternionic groups. 

(2) 

Now it is still not obvious that this group is an invariance group of C D , where 
the superscript ^ indicates the presence of two fermionic fields (with the same mass) 
since the 7^ matrices are complex, not real. However, it is readily verified that an 
equivalent set of 7^ matrices (with the same commutation relations) are given if the 

2 We recall the well known fact that Cd becomes identically null, if the field equations are applied 
to tp and 
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z-factors in 7 M are substituted by 1 | i. Whence the imaginary units appear to the 
right of ip and the elements given by eq. ( p5[ ) commute with 7 M and hence cancel 
within Cd- One may also derive the similarity transformation which performs the 
above change of representation and thus derive the corresponding representation^ of 
{7(1, q) for our original (conventional) choice of complex 7 M . 

We note that our {7(1, q) group consists of quaternionic "phase" factors multiplied 
by the unit 4x4 matrix^. We may ask if there are other 4x4 (essentially real) matrices 
which commute with the four dimensional Dirac gamma matrices. However there is 
no other matrix that commutes with all four gamma matrices so the maximum global 
invariance group is indeed {7(1, q) | {7(1, c) to be identified with the (complex) 
Glashow group SU{2) x {7(1). 

5 Conclusions 

We begin our conclusions from the end results of the last section. We wish to recall 
that not every global invariance group is automatically gauged into the corresponding 
local group in gauge theory. For example, the so-called complex bosonic Lagrangian 
may a priori (before eventual gauging) represent two neutral equal mass particles or 
a complex pair of oppositely charged particle-antiparticle. Thus even in the case of 
a Dirac Lagrangian with complex (4 x 4) Dirac matrices we do not have necessarily 
a doublet under {7(1, q). So that we have not derived the gauge group of the Salam- 
Weinberg model. The pair may consist of two singlets. However it was not even 
obvious a priori that a global invariance group {7(1, q) isomorphic at the generator 
level with SU(2, c) existed. We have thus shown that even with quaternionic fields, 
it is possible to impose a Glashow group invariance and that this occurs by merely 
adopting reducible gamma matrices. 

Our viewpoint is that SU(2, c) invariance in particle physics is really indication 
of invariance under the simplest of all unitary quaternionic groups. We also have 
a complex invariance group {7(1, c) but this is justified by the complex projection 
of the Lagrangian density. This complex projection was not imposed to obtain the 
{7(1, c) factor group but in order to obtain the desired quaternionic Dirac equation. 
The fact that this group exists in nature as the weak-hypercharge group is an un- 
doubted success of this model. Further, the automatic appearance of this complex 
unitary group is expected whatever the left acting (quaternionic) unitary group is. 
This strongly suggests that in the search for grand unified theories one should con- 
sider preferentially a product group of the type G \ {7(1, c) with G a quaternionic 
unitary group, conditioned to commute with the Dirac gamma matrices in the chosen 
representation. 

3 We observe that while our {7(1, q) group acts from the left on tp and {7(1, c) acts on the right, 
this has nothing to do with the helicity indices L and R of the Weinberg-Salam theory. 

4 An additional right acting complex phase {7(1, c) is also allowed, but this is equivalent to that 
already described in the previous section. 
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Let us recall the other results of this paper. We have discussed the application of 
the variational principle to Lagrangians with possibly quaternionic fields. We noted 
en passant that even within standard field theory the limitation of variations in the 
fields to complex variations is an implicit assumption, since the variations themselves 
have no physical content. Thus quaternionic variations have always existed even if 
not traditionally applied. Of course this observation can be generalized to even more 
exotic variations, e.g. supersymmetric variations containing Grassman terms. 

Finally, whereas in this work the need of a complex projection of the Diarc La- 
grangian density is demonstrated, for the scalar field Lagrangian it is an assumption 
with interesting consequences ||. We suggest that, for reasons which still elude us, all 
terms in the quaternionic Lagrangians including interaction terms must be complex 
projected. 
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